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We perform a study of the gravitating electrostatic spherically symmetric (G-ESS) solutions of 
Einstein field equations minimally coupled to generalized non-linear abelian gauge models in three 
space dimensions. These models are defined by lagrangian densities which are general functions of 
the gauge field invariants, restricted by some physical conditions of admissibility. They include the 
class of non-linear electrodynamics supporting ESS non-topological soliton solutions in absence of 
gravity. We establish that the quafitative structure of the G-ESS solutions of admissible models 
is fully characterized by the asymptotic and central-field behaviours of their ESS solutions in flat 
space (or, equivalently, by the behaviour of the lagrangian densities in vacuum and on the point 
of the boundary of their domain of definition, where the second gauge invariant vanishes). The 
structure of these G-ESS configurations for admissible models supporting divergent-energy ESS 
solutions in flat space is qualitatively the same as in the Reissner-Nordstrom case. In contrast, the 
G-ESS configurations of the models supporting finite-energy ESS solutions in flat space exhibit new 
qualitative features, which are discussed in terms of the ADM mass, the charge and the soliton 
energy. Most of the results concerning well known models, such as the electrodynamics of Maxwell, 
Born-Infeld and the Euler-Heisenberg effective lagrangian of QED, minimally coupled to gravitation, 
are shown to be corollaries of general statements of this analysis. 

PACS numbers: 04.40.-b, 04.70.Bw, 05.45.Yv, ll.lO.Lm 



I. INTRODUCTION 



For more than two decades a great deal of attention 
has been focused on the study of gravitating field con- 
figurations. This renewed interest is partially motivated 
by the discovery that generalized field theories, such as 
the Born-Infeld (BI) lagrangian (introduced in the 30's 
to remove the divergence of the electron self-energy in 
classical electrodynamics and some non-abelian BI- 
like versions arise, together with the gravitational field, 
in the low-energy limit of string theory [2] . But the in- 
terest in solutions of the Einstein equations coupled to 
different kinds of fields has other motivations. Indeed, in 
a four-dimensional flat space, several non-existence the- 
orems [3] restrict drastically the class of field theories 
supporting soliton solutions. However, through the cou- 
pling to gravity this obstruction can be removed and thus 
gravitating particle-like solutions can be found in theo- 
ries that do not support soliton solutions in flat space, 
such as the pure Yang-Mills one On the other hand, 
theories supporting soliton solutions in flat space, as the 
abelian and non-abelian BI models, have been extended 
to curved space leading to black hole-like solutions (see 
for a review). It is also worth mentioning the so- 
called soliton stars, gravitating coherent quantum field 
states with the features of non-topological solitons [)§] . 

In the context of electromagnetic field theories, gener- 
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afizations of the Reissner-Nordstrom (RN) solution of the 
Einstein-Maxwell field equations have been developed for 
some particular examples of non-linear electrodynamics 
(NED) coupled to gravity. The first finding in this field 
dates back to the work of Hoffmann [7]. In the eighties, 
aside from the electrically charged black hole solutions 
in the Einstein-BI system [8\ , similar configurations were 
found in more general Bl-type electrodynamics min- 
imally coupled to gravity. These studies were later ex- 
tended by the inclusion of a dilaton field [iO|jthe anal- 
ysis of black holes in Anti-de Sitter spaces 11| and its 
generalization to higher-dimensional cases |12| . Besides 
Bl-type fields the coupling to gravity of other NED la- 
grangians has been also considered. Let us mention the 
gravitating electrostatic, monopole and dyon solutions of 
the Euler-Heisenberg (EH) effective lagrangian of QED 
[Ts! ]. the logarithmic lagrangian of Ref.[14] and the class 
of models obtained as powers of the Maxwell lagrangian, 
studied in Ref . . 

In the search for regular (black hole or not) solutions, a 
theorem [l^ forbids their existence for electrostatic fields 
in NED models having the Maxwell weak field limit. But 
this theorem can be circumvented for purely magnetic so- 
lutions [l3| , or through a coupling between different fields 
[isj . Moreover, some gravitating NED models were re- 
ported to lead to electrically charged non-singular black 
hole solutions lo*], but they are based on unphysical 
(multivalued) lagrangian density functions (20| . 

In this work we perform a general study of G-ESS solu- 
tions for a large class of NED models minimally coupled 
to gravity. The lagrangian densities for these abelian 
gauge fields are defined as arbitrary functions of the two 
standard gauge invariants and include most of the afore- 
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mentioned models. This class is restricted by some phys- 
ically reasonable "admissibility" conditions, such as the 
regularity and uniqueness of the lagrangian densities, the 
positivity of the energy, the parity invariance and the 
asymptotic vanishing of the ESS fields (faster than 1/r) . 
These admissible models can be exhaustively classified 
into two sets. The first one includes the models with 
ESS solutions in fiat space which are energy-divergent, 
owing to their central-field behaviours (the Maxwell the- 
ory being the simplest example). The second one includes 
the models supporting finite-energy ESS solutions in fiat 
space, which were extensively analyzed in Ref.[2l|. The 
structure of the G-ESS solutions of the first set is shown 
to be the same as that of the RN solutions of the Einstein- 
Maxwell field equations. The "soliton-supporting" NED 
models, when coupled to gravity, lead to new qualita- 
tive features of the G-ESS solutions (single-horizon black 
holes, "black points" , etc) which we analyze and classify 
in terms of the charge, ADM mass and flat-space soliton 
energy. We shall show that the central-field behaviours 
of the corresponding ESS solutions in flat space which, as 
already mentioned, allowed the determination and clas- 
sification of the full set of this kind of admissible models 
[2l| , allow also the complete characterization of the corre- 
sponding gravitational configurations. This central-field 
behaviour is determined by the specification of the form 
of the admissible lagrangian densities on the boundary 
of their domain of definition for vanishing values of the 
second gauge invariant. Consequently, the gravitational 
structure of these solutions is fully characterized by ad- 
missibility and this boundary behaviour, regardless of the 
explicit expressions of the lagrangian densities elsewhere. 
The G-ESS solutions for some models studied in the liter- 
ature (as the gravitating BI 8] and EH 13]) are exposed 
as representative examples of the different classes intro- 
duced here, and their gravitational features are shown to 
be immediate consequences of the general results derived 
in the present analysis. 



the fine {X > 0,Y = 0)f]n a.nd dip/dX to be strictly 
positive there. We also require <f{X, Y) to be symmetric 
in the second argument {'fiX, Y) = (p{X, —Y)), in order 
to implement parity invariance. Moreover, the require- 
ment of the positive definite character of the energy is 
essential in the present context and leads to the minimal 
necessary and sufficient condition on the energy density 
p, (obtained from the symmetric energy-momentum ten- 
sor of Eq. Q below) 

p> (v/x^Ti^ + x)|| + r|^-^(x,y)>o, (1) 

to be satisfied in the entire domain of definition fl. 

These models are minimally coupled to gravity through 
the action 
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where, as usual, R is the scalar of curvature and g is 
the determinant of the metric tensor g^^. The associated 
Euler field equations, together with the Bianchi identities 
for the electromagnetic field, take the form 



(3) 



where ipx — §^ and ipy — On the other hand, 

the gauge-invariant symmetric energy-momentum tensor 
reads 



—2 SSned _ 
2 (ifi^F^^F^ - VyF^o^F^) ~ g^^MX.Y). (4) 



II. MODELS AND FIELD EQUATIONS 



For ESS solutions we have 



Let us precise our conventions and conditions. In 
(3-1-1) dimensions the two quadratic invariants of the 
abelian field, built from the Maxwell field strength ten- 
sor Ff,^ = df,A^ - d^Af, and its dual F* = ^e^^^/jF"'', 



are defined as X 



- and Y 



— i^^j^F*^" — 2E ■ B, where the electric and magnetic 
fields are F* = — F"* and B^ = —^e^^^Fjk, respectively. 
The dynamics of these fields is governed by lagrangian 
density functions ip{X^ Y) of these invariants. For phys- 
ical admissibility [21] these functions must be single- 
branched, continuous and derivable on their domains of 
definition (fi S 3?^) of the X — Y plane, which are as- 
sumed to be open, connected and including the vacuum 
{{X = 0,y = 0) e n). The regularity (for r 7^ 0) of the 
ESS solutions requires also ip{X, Y) to be of class on 



F : 



F(r)- 
r 



B = 0. 



(5) 



In this case the energy-momentum tensor (j4]) satisfies 
Tq — Tl and, as a direct consequence, the metric can be 
cast into the Schwarzschild-like form 



where dVi^ = d9^ + sin 



dr^ 



r^dn^. 



(6) 



In this coordinate system 
the non- vanishing components of the symmetric energy- 
momentum tensor for ESS fields read 
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from which the Einstein equations are immediately writ- 
ten. Under the conditions of admissibility stated above 
(in particular, Eq.([T|)) it can be easily shown that the 
Weak Energy Condition T^jy^''^"^ > (where is a time- 
like vector) automatically holds. This condition ensures 
that a time-like observer measures a non-negative energy 
density of the field at any point. 

The generalized Maxwell equations ([3]) for these ESS 
fields take now the same form, in the Schwarzschild co- 
ordinate system, as the corresponding equations in the 
flat space problem (in spherical coordinates) and have 
the same first-integral, which reads 



(8) 



where q is an integration constant, identified as the elec- 
tric charge [2l| (without loss of generality, owing to the 
positivity of (fxO'nY — 0, which is a consequence of 
Eq.(IT]), we consider the case q > 0, E{r) > only). The 
invariant X is now given hy X — E'^ (note that the sec- 
ond set of equations in ^ is identically satisfied for these 
fields). Let us emphasize the fact that the first- integral 
(|8]) does not depend explicitly on the metric function A(r). 
It determines implicitly the form of the field E{r), once 
the expression of the lagrangian function ^p{X^ Y) is spec- 
ified. For a given lagrangian this form is the same, in the 
Schwarzschild coordinate system, as in the absence of 
gravitation in spherical coordinates. Moreover, the ex- 
pression of TQ{r) (as well as the other components of the 
energy-momentum tensor) in ([T]) is also the same as in 
the flat space case. Consequently, the results of the anal- 
ysis of Ref. 2l| on generalized NEDs can be immediately 
translated to this gravitating problem. Let us recall some 
main results of this reference which are pertinent for the 
present study. We have characterized the lagrangians of 
the NEDs by the behaviours of their ESS solutions at 
the origin and at inflnity, in order to establish the con- 
ditions for the convergence of the integral of energy. Let 
us assume in both limits a behaviour of the form 



p/2 

E{r) ^ ( ^ 1 r- 



(9) 



(with p 7^ 0), /3 being a characteristic constant for each 
model, given by the limits 



/3^1im(X^|| 



(10) 



calculated as X -> or as X ^ oo for the asymptotic 
or central-fleld behaviours, respectively. We imposed 
p<— lasr— >oo (in order for the integral of energy 
to converge as r — > oo) and we distinguished three cases, 
corresponding asymptotically to a slower than coulom- 
bian damping (—2 < p < —1, case Bl ), a coulombian 
damping (p = —2, case B2 ) and a faster than coulom- 
bian damping {p < —2, case B3 ) of the fleld. As r — )■ 



we considered the models with —1 < p < (the only 
admissible ones compatible with the convergence of the 
integral of energy at the center) and we distinguished two 
cases. The case Al (— 1 < p < 0), which corresponds to 
ESS fields divergent (but integrable) at the center, be- 
having as in Eq.®, and the case A2 {p = 0), which cor- 
responds to a constant value of the field at the center, 
where it behaves as 



E{r 0) 



br" 



(11) 



a and a being positive constants characteristic of the 
model (for example, for the BI model, a — A and the 
maximum field strength is E{0) = a). The constant b 
is related to the charge of the particular ESS solution 
through 



a = lim 



a{a 



.2/a 



dip 
dX 



= qb^/" 



(12) 



a being a characteristic parameter of the model. 

The behaviours of the lagrangian densities ipiX, Y — 
0) corresponding to the asymptotic behaviours of the 
finite-energy ESS solutions take the generic form (see 
Fig.l) 



^{X,Y = 0)^^X\ 
7 



around the vacuum {X = E'^{r 



(13) 

oo) ^ rP ^ 0), 
where 7 = ^i^- ^'-'^ integral of energy to converge 
asymptotically this parameter must range in the interval 
1/2 < 7 < 3/2, with 1 < 7 < 3/2 in case Bl , 7 = 1 in 
case B2 and l/2<7<lin case B3 . 

The same expression (|13p with 7 > 3/2 describes 
the behaviour of the lagrangian density for large X in 
case Al . This function diverges on the boundary of 
at the point Y — and X = E'^{r ^ 0) rP ^ 00, 
exhibiting there a vertical parabolic branch. In case A2 
the boundary of at F = is the point {X = E'^{r = 
0) — ,Y = 0) and the lagrangian density behaves there 
(for (7 ^ 2) as 



ip{X) 



2aa 



ia-VX)—+ip{a\Y^O), 



(14) 



exhibiting a vertical asymptote in this point for cr < 2, or 
taking a finite value there with infinite slope for a > 2. 
For a — 2 the lagrangian density behaves as 



(p{X) —a ln(a 



(15) 



and exhibits a vertical asymptote ai X = a^. 

In putting aside these conditions we are lead to six 
patterns of models for the class of admissible NEDs sup- 
porting (flat-space) finite-energy ESS solutions. Never- 
theless, here we are also interested in the gravitating ver- 
sions of NEDs supporting energy-divergent (due to the 
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central field behaviour) ESS solutions in flat space, as 
the RN solution of the Einstein-Maxwell field equations. 
Consequently, in the sequel we shall consider also cases 
where the ESS field solutions diverge at the center faster 
than l/r {p < —1 in Eq.® for r 0, corresponding to 
1/2 < 7 < 3/2 in Eq.dTSl) for large X). 

(p{X,Y=0) 
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FIG. 1: Behaviours of the lagrangian densities of admissible 
NEDs for y = 0. Around the vacuum (X — 0) they determine 
the asymptotic behaviours of the ESS solutions, necessary for 
the integral of energy to converge at infinity (cases Bl, B2 and 
B3). The central behaviour of the ESS solutions character- 
izes the structure of the corresponding G-ESS configurations. 
This behaviour is governed by the form of the lagrangian den- 
sities on the boundary of their domain of definition. The 
cases Al and A2 correspond to ESS solutions for which the 
integral of energy converges at the center. In case A2 the 
boundary is reached at finite X and the lagrangian densities 
exhibit there a vertical asymptote (if cr < 2) or an absolute 
maximum with divergent slope (if a > 2). In case Al the 
boundary is at infinity and the lagrangian densities diverge 
asymptotically as (p{X — > oo) --^ X'' with 7 > 3/2, exhibiting 
a vertical parabolic branch there. Lagrangian densities with 
similar asymptotic form and 1/2 < 7 < 3/2 {p < —1), lead- 
ing to ESS fields with a divergent integral of energy (at the 
center), have also been drawn. 

Note that the field equations can now be deduced from 
the Bianchi identity V^G'^'^ = and, consequently, ([8]) is 
also a first-integral for the Einstein equations. By com- 
bining Einstein equations and the first-integral we are 
lead to the system (we take G = 1) 



— (rA(r) - r) 



(rA(r)) ~ — IGTrrTf ~ IQnrip, 



(16) 



whose compatibility can be easily established. Thus the 
solution of the G-ESS problem reduces to solve the first 
of Eqs. ([TB| . when E{r) is determined from Eq.([5]), once 
the form of the lagrangian function is specified. The in- 
tegration of leads to the expression 



A(r,g) = 1 



2M , 2£ex(r,g) 



(17) 



where M = i \imr^oo[r{l - A(r))] is an integration con- 
stant, identified as the ADM mass, and 



POO 

£ex(r,(z) =4^ / R^T^{R,q)dR, 

J r 



(18) 



is the energy in flat space outside the sphere of radius 
r (by convenience we shall call "integrals of energy" the 
expressions similar to (dH)). The integrand of ([T5)) in 
terms of the held E{r,q) can be obtained from Eqs.([7]) 
and ^ and reads 



(19) 



2q 



E{r,q) 



dE{R,q)/dR 
R^ 



dR 



whose asymptotic behaviour can be shown to be the same 
of E{r, q). Then the integral deflning SexiT, q) is conver- 
gent for r > 0. The solution PT)) contains the charge q 
and the mass M as arbitrary constants. The dependence 
of the integral of energy on the charge is given by 



'Xr,q) 



- d^l'^e 



(20) 



Another important result is that, for any admissible 
model, Eea; ij , q) is a monotonically decreasing and 
concave function of r. This can be proven by dou- 
ble derivation of ([T8| . taking into account (fT9|) and the 
monotonically decreasing character of E(r), which can 
be established from ([5]) and the admissibility conditions 
(see Ref.[2l[). Another useful expression for the exterior 
integral of energy in terms of the electrostatic potential 
can be obtained from Eqs. ([T8|) . (fT9|) and ([8|) and reads 



'-^Mr,q)^'-^r^T^ir,q). 



(21) 



All these results allow the complete characterization 
of the structure of the G-ESS solutions of any physically 
meaningful gravitating NED of the form We shall 
begin by the analysis of the family of those models lead- 
ing to energy-divergent ESS solutions in flat space. 



MODELS WITH ENERGY-DIVERGENT 
ESS SOLUTIONS IN FLAT SPACE 



III. 



Let us return to Eg. pTj) and look for solutions rh of 
A(r) = 0, which give the radii of the horizons present in 
the configuration. In obtaining the solutions r;i 7^ we 
can solve the equivalent equation 
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M 



(22) 



which simphfies greatly the discussion and will give the 
inass-horizon-radius relation for the black holes. Indeed, 
the solutions of (|22|) are the intersection points between 
the curves y — e^^ (r, q) and the beam of straight lines 
y ~ M—r/2 in the r—y plane (see Fig. 2). The monotonic 
character and the concavity oiy — Sex (f, q) for any value 
of q lead to three possibilities: 

(1) There is a family of extreme black hole solutions 
(with parameter q) which correspond to the points of the 
exterior energy curves with slope —1/2. The correspond- 
ing radii of the horizons are the solutions rhextril) of the 
equations 



8^r2TO(r,<7) = l, 



(23) 



and the associated values of the ADM masses, Mextr{q), 
are given by replacing rhextril) into the general mass- 
horizon radius relation (|22p . Using Eqs. ipTj) and we 
are lead, in terms of the electrostatic potential, to the 
extreme black hole mass formula 



Mextriq) = 



rhextriq) , 167rgf 



which is valid for G-ESS fields vanishing asymptotically 
faster than 1/r. As easily shown from Eq. (|23| , the scaling 
law ([20| and the shape of Sex {r) , as q increases from 
to 00, rhextr{<l) and M^xtril) get increased in the same 
range. 




(2) For configurations with values of M < Mextril) 
there are no horizons and the solutions have a naked 
singularity at the center. 

(3) For configurations with values of M > Mextr{q) 
there are always two horizons (event and Cauchy) for 
each configuration. The corresponding radii are the so- 
lutions of ([22|) . In both (2) and (3) cases, the limit 
M — )• Mextriq) is singular (see Ref.|22|V 

Let us stress that all these results on the gravitational 
properties of the ESS solutions at finite r are indepen- 
dent of their asymptotic behaviours. They are similar to 
those of the RN solution, corresponding to the Maxwell 
lagrangian density 



ip{X,Y)^aX, 



(25) 



(a being a constant) which characterizes the pattern for 
the structure of the gravitational field of the solutions of 
this family. In this case Eqs. ([22|) and ((23|) can be solved 
explicitly and, aside from Eqs.© and (EL lead to the 
well known results (see, for example, Ref.[23j) 



E{r,q) 
rh{q) 



r^T^{r,q) = ; eex{r,q) = 



ar 



+ ^Mr.exUqU), (24) 




ar 
(26) 



IV. MODELS WITH FINITE-ENERGY ESS 
SOLUTIONS IN FLAT SPACE 

In this case the integral of energy converges at the 
center (r — > 0) and we can define both the exterior and 
the interior integrals of energy for a given sphere 

eUr,q) = 4^ / R^T^{R,q)dR 
Jo 

poo 

eex{r,q) = 4^ / R^T°{R,q)dR, (27) 

J r 

which are related through 

siq) = £in{oo, q) = £ex{0, (?) = £in{r, q) + £ex{r, q), (28) 

where e{q) is the total fiat-space energy of the ESS solu- 
tion. The useful formula 



FIG. 2: Exterior integral of energy esxXr) for NEDs with 
(flat-space) energy-divergent ESS solutions. The straight 
lines correspond to different values of the ADM mass M. 
The cut points give the horizons of the black hole solutions: 
(1) tangency point corresponding to single-horizon extreme 
black holes (M = Mextriq))', (2) naked singularities (M < 
Mextr{q)); (3) black holes with two horizons (M > Mextr{q))- 



s{q) 



IGtt 



,3/2 



-1/2 



dy, 



(29) 

allows the explicit calculation of this total energy, once 
the expression of the lagrangian density is known [2l|. 
Using Eq. ipn]) we obtain the scaling law 
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Siq) = q'/^iq = 1) 



(30) 



where e(g = 1) is the flat-space energy of the unit charge 
solution and is a universal constant for a given model. In 
terms of the interior integral of energy, equation (jl7p can 
be alternatively rewritten as 



limit case, which steps from +00 to —00 on both sides of 
the vanishing-radius inner horizon. In both (Al-2) and 
(Al-3) cases the limit M — > Mextr{q) is singular. 

(Al-4) When M = e{q) we have (see Eqs.jni) and 

m) 



A(r,g) = 1 



2ein{r,q) 



-00 



as 



r ^ 0, (33) 



A(r,q,M) = l 



2{M-e{q)) 2e„,ir,q) 



(31) 



The presence of this finite energy, as a parameter to 
which the ADM mass can be compared, leads to the exis- 
tence of special solutions for which the total energy of the 
configuration, as defined by the ADM mass, equals the 
electrostatic energy in absence of gravitation. We shall 
call these solutions "critical" . 

The function Sexi^iq), for fixed q, is again monotoni- 
cally decreasing and concave, whereas Sin(r,q) is mono- 
tonically increasing and convex. Returning to Eq. (jl7p we 
see that the roots of X{r,q,M) = 0, giving the radii of 
the horizons for fixed q and M, can be obtained (when 
non- vanishing) by solving the same equation (|22p . But 
now the function eex(r, q) reaches the finite value e{q) at 
r = 0. Let us analyze separately the cases Al and A2. 

Case Al : Near the origin the field behaves as 

E{r,q) ^ r]{q)rP, where — 1 < p < and r]{q) — (^/^^ 

(see Eq.®). From Eg. lfTn)) we obtain for eexir^q) the 
behaviour 



c{r,q) ^ e{q) 



16Trqri{q) 
{2-p){l+p)' 



(32) 



around the center, where its slope diverges. The study 
of the intersection points of the curve y — Sex (f, q) (for 
fixed q) with the beam of straight lines y — M — r/2 
gives the masses and radii of the horizons and leads to 
the following classification of the solutions (see Fig. 3): 

(Al-1) There is again a family of extreme black hole 
solutions (with parameter q) whose radii rhextr{q) and 
ADM masses Mextr(q) are given by the same equations 
([^ and (1231) as in the divergent-energy case of section 
III. These are increasing functions of the charge, ranging 
from to 00 as q increases in the same interval, as can 
be shown from the scaling law (|20l) for the integral of 
energy. 

(Al-2) Again, for configurations with values of M < 
Mextriq) there are no horizons and the solutions have a 
naked singularity at the center. 

(Al-3) For Mextriq) < M < e{q) the solutions exhibit 
two horizons. As M — e{q) — >■ 0~ the radii of the inner 
horizons of the sequence of solutions parameterized by 
M goes to zero and the radii of the outer horizons ap- 
proach a regular limit value rhcritiq), given by the non- 
vanishing solution of the equation r = 2£i„(r, g). There 
is a singularity of metric function A at the center in this 



and only the outer horizon of radius ThcritiQ) survives 
(see Fig. 4). These solutions coincide, for r > 0, with 
the limit solutions of the sequence of the previous case 
(Al-3), but differ by the nature of the singularity at the 
center. 

(Al-5) For M > e{q) there is a unique horizon whose 
radius, obtained from Ea.p2|). increases with M from 
fh — fhcrit{q) when M — s{q), becoming ^ 2M for 
large M. As M — e(g) — >■ 0+ (for fixed q) the sequence 
converges continuously to the critical solution of the case 
(Al-4) for any r. 

As an example of this family let us consider the effec- 
tive action for QED of Euler and Heisenberg j2^] , whose 
lagrangian density takes the form 



(34) 



where /i is a positive constant. This model satisfies the 
admissibility conditions and is a particular case of a fam- 
ily of polynomial lagrangian densities, all belonging to 
this class Al, whose ESS solutions in fiat space were ex- 
tensively analyzed in Ref.[2l|. The ESS solutions in this 
case are obtained from Eq.®, which now takes the form 



2^£;3(r,q) + ii?(r,g) = \, 



(35) 



and can be solved explicitly through the Tartaglia for- 
mula, leading to 



E{r,q) 







1/3 
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1/3 



(36) 



where A 



-I- if^ > 0, u 



rrk- and V 



Near the 



center these fields behave as E{r 0,q) ^ i^Y^'^''" ^^'^ 
(case Al) and are asymptotically coulombian: E{r 
00, q) ^ ^ (case B2). These behaviours can be also 
immediately deduced by the simple inspection of the la- 
grangian density ^ and Eq.®. The central field be- 
haviour {p = —2/3), together with the admissibility con- 
ditions endorse the decreasing and concave character of 
the exterior integral of energy. Consequently, the state- 
ments Al-1 to Al-5 must hold for the gravitating version 
of this model. Let us calculate some characteristic mag- 
nitudes in order to obtain the quantitative description of 
the G-ESS configurations. 
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The energy density obtained from Eqs.Q and (p4|) can 
be written, in terms of the field, as 



(37) 



and its central and asymptotic behaviours are 



TO(r^0,q) 



3<Z 



,4/3 



(16a.)V 



(38) 

confirming that the slope of the integral of energy 
£ex{r,q) diverges as r — (see Fig 3). 

Using Eq. (f29| the total flat-space energy associated to 
these solutions takes the form 



mnq' 



,3/2 foo 



dy 



/o ^y(l/2 + 2/iy2) 3^*1/4 A' 

(39) 

where B{x,y) — t^^^{l — ty~^dt is the Euler inte- 
gral of first kind. The extreme black hole solutions must 
satisfy Eg. ([25)1 . which now reads 



leading to the relation 



2 



(40) 



function of q in the same interval, in agreement with the 
statement Al-1. Moreover, Eqs. (H51) and (H^ define ex- 
plicitly the function q{rhextr), giving the exact relation 
between the charge and the horizon radius of the extreme 
black holes. The associated mass Mextr{q) can now be 
obtained from Eq. ([M|) (after integration of the field p6l) . 
leading to the form of the potential Ao{r, q)) and, as eas- 
ily seen, it is also an increasing function of q or VhextriQ), 
ranging from to oo. It is now easy to verify that so- 
lutions with M < Mextriq), Mextriq) < M < e{q) and 
M > £{q) lead to naked singularities, two-horizons black 
holes and single-horizon black holes, respectively, the pre- 
vious equations allowing the quantitative determination 
of the parameters of these configurations. 



y 



£ (q 




between the radii of the horizons of the extreme black 
holes and the values of the electrostatic field there. Let 
us introduce, by convenience, the auxiliary variable 

X=Jl+ 2^^,, -l> (42) 

which is a monotonically decreasing and concave function 
of Thextr , ranging from oo to as r/ieatr ranges from to 
oo. Another independent relation between E^xtriq) and 
X is given by the first integral (j35p . and the elimination 
of Eextriq) between both equations leads to 



FIG. 3: Exterior integral of energy eex{r,q) for the two dif- 
ferent cases of NEDs with (flat-space) finite-energy ESS so- 
lutions. The bottom curve corresponds to models with ESS 
fields divergent (but integrable) at the center (case Al). In 
this case the slope of the curve at the origin diverges. The 
remaining curves are representative of ESS solutions of NEDs 
leading to a finite value of the fields at the center (A2 cases). 
They correspond to the slopes at the center associated to dif- 
ferent values of the electric charge {e'exi''' = 0) = —Siraq ^ 
1/2, see Eg. (144^ 1. The cut points with the straight (dashed) 
lines y = (M— r/2) give the horizons of the different black hole 
solutions (see the discussion in the text). For the sake of sim- 
plicity curves corresponding to different values of the charge 
q have been superimposed on the same figure, in such a way 
that all of them begin at the same point e{q). Consequently, 
each curve should be viewed as if it were represented in a dif- 
ferent figure. Note that the slope of the beam of straight lines 
is not affected by this trick. 



2„2 



Un-'q 



-X 



1 



X- 



(43) 



By drawing the functions of x defined by the both sides 
of this equation (a beam of straight lines with slope pro- 
portional to q^, cutting once a squared hyperbola) it is 
obvious that there is a unique solution xiq) > 0, which 
is a decreasing function of q, ranging from oo to as 
q ranges from to oo. Hence rhextriq) is an increasing 



Case A2 : As above mentioned, this case corresponds 
to p = near the origin and the ESS field solutions be- 
have there as in Eg. pTj) . Before using Eq. ip^ in order 
to analyze the structure of the G-ESS solutions we must 
study the behaviour of eex{i',q) around r — 0. Using 
Egs. p^ . (P7)) and pT|) we obtain the expression 



dr 



(44) 
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for the slope of £ex{T^q) at the center. Moreover, the 
behaviour of the gravitational potential around the center 
can be obtained from Eqs. (f3T|) . (fT9|) and (fTTj) and reads 



A(r) ~ 1— 167raq^ 
for cr 7^ 2 and 



{a + l){2-a) 



(45) 



AW 



1 — 167ra(7 



2{M - siq)) 



+ 



^r2(l-21n(r))+Ar2, 



(46) 



for a = 2, A being an integration constant. When a > 2 
this constant is given by A = 87rgy(a ,o) ^^^^ Ea.([T4|V 
Otherwise the value of A is not relevant because the asso- 
ciated term is not dominant. Consequently, we must now 
analyze several possibilities corresponding to solutions 
with different values of q, which we shall denote as A2a 
(corresponding to g < l/167ra), A2b (for q > l/lGira) 
and A2c (for q = l/167ra). Let us discuss the behaviour 
of the configurations corresponding to different values of 
M for each one of these cases. 

• (A2a) If Wnaq < 1 we have three classes of solutions 
(see Fig.3): 

(A2a-1) M < e{q): The solutions have no horizons 
but only naked singularities at the center. As M—e{q) -> 
0~ the solutions, parameterized by M, converge towards 
a limit field with a discontinuity of A(r) at the center, 
which steps from +oo to the finite value 1 — lQT:aq > 0. 
For r > 0, A(r) converges towards the form of the critical 
case (A2a-2) below. 

(A2a-2) M = e{q): In this case A(0) = 1 - 16™^ > 
is finite and A(r) remains regular and positive every- 
where. The behaviour of the derivative of the metric near 
the center can be calculated from the expression (|45l) and 
reads (for (7^2) 



dX 

dr 



32TT(7bq 



(a + l)(2-a) 



2rA, 



(47) 



which diverges for cr < 1, vanishes for a > I and takes a 
positive finite value for cr = 1 (see Fig. 4). Equation (j46|) 
shows that this derivative also vanishes for a = 2. These 
are naked singularities whose structure at r = differs 
from that of the limit field defined in (A2a-1), but both 
fields coincide for r > (see Fig. 4). 

(A2a-3) M > e{q): The solutions have a unique and 
non-degenerate horizon. The radius of the horizon rh 
is related to the ADM mass through Eq.(l22]). At the 
center A(0) = -oo. As M - e{q) ^ 0+, r/i ^ and 
the solutions, parameterized by M, converge towards a 
kind of ^^black poinf\ with a discontinuity of A(r) at the 
center, which steps from — oo to zero on the internal side 
of the horizon {r^), and from zero to 1 — IGnaq on the 



external side (r^ ) . For r > this limit solution is regular 
and coincides with the one of the critical case (A2a-2). 

• (A2b) When IGnaq > 1 we have five classes of solu- 
tions (see Fig.3 and the small frame of Fig. 4): 

(A2b-1) Now, as in the case Al, there is a family of 
extreme black hole solutions (with parameter q) whose 
radii rhextr{q) and ADM masses Mg^triq) are the so- 
lutions of Eqs. (|23l) and (|24|) (see curve B in the small 
frame of Fig. 4). The analysis of these equations, using 
the scaling law (|20l) for the integral of energy, shows that 
Thextriq) and Mextr{q) are monotonically increasing func- 
tions of g, ranging from to oo and from ^(y^^) to oo, 
respectively, as q ranges from q 



167ra 



to OO. 



(A2b-2) For M < Mextr{q)'- There are no horizons 
and the solutions have a naked singularity at the center, 
analogue to the one of case (Al-2). As in the cases of 
RN and (Al-2), the limit M — Mextr{q) is singular (see 
curve A in Fig. 4). 

(A2b-3) For M^^tM) < M < e{q): As in case (Al- 

3) the solutions have two horizons. As M — e{q) 0^ 
the radii of the inner horizons vanish and the sequence 
of solutions converges to black holes with an external 
event horizon and an infinite jump of the metric at the 
center, which passes from -|-oo to the finite value A(0"'") = 
1 — 16naq < on both sides of the (vanishing-radius) 
inner horizon (see curves C and D in Fig. 4). As in the 
RN case, the limit M — Mextrio) 0+ is singular. 

(A2b-4) When M ^ £{q): As in case (A2a-2) A(0) = 
1 — IQ-Kaq < is finite (but now negative), and A(r) 
is regular everywhere and vanishes at some r — rhcrit, 
where there is a unique horizon. The derivative of the 
metric near the center takes the same form (|T7l) as in the 
(A2a-2) case, with the same behaviour as a function of cr. 
These are now black holes which coincide, for r > 0, with 
the limit cases M—e{q) — 0~ of the sequence of solutions 
(A2b-3) (it coincides also with the limit of the sequence 
(A2b-5) below, as M - e{q) -J> 0+) but the structure of 
the singularities at the center is different in each case (see 
curves D and E in Fig. 4). 

(A2b-5) For M > e{q): There is now a unique hori- 
zon with a radius ranging from = Vhcrit to oo as 
M increases. The sequence of solutions converges as 
M — e[q) — >■ 0+, for r > 0, to the critical case (A2b- 

4) , but now A(r) steps from — oo to the negative value 
A(0) = 1 — l&iraq < at the center (see curves E and F 
in the small frame of Fig. 4) . 

• (A2c) When l&Traq = 1 the charge is fixed, Eq. p5)) 
is satisfied at r = and Eqs. (|^51) and (|46p get simplified. 
We can distinguish three cases (see Fig.3): 

(A2c-1) For M ~ £{q)- In this case A(r) is regular 
everywhere and vanishes at r = 0, where the condition for 
the radius of extreme black holes ((23|) holds. There are 
three behaviours for the slope of A(r) at r = (depending 
on the values of cr = 1) given by Ea. (|47)) . These critical 
solutions are extreme black points (see Fig. 4). 

(A2c-2) For M < e{q): We have naked singularities 
with the same behaviour as in case (A2a-1) but in the 
limit M — e{q) — > 0~ an extreme black point appears 
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with a discontinuity of A(r), which steps from +00 to 
zero at the center. This limit is singular. 

(A2c-3) For M > e(q): The behaviour is similar as 
in (A2a-3) and we have a unique horizon. At the center 
A(0) = —00, and the solutions converge to black points 
as M — e{q) — > 0+. Once more the limit of the sequences 
coincides with the solutions of the critical case (A2c-2) 
for r > 0, but they have point-like discontinuities at the 
center. 



(A2b) 



(Al) 



a<l 



a>l 



1 \ \M<£(q) 


\ 






'FT 

FAhe 


k ^ 'a 

/ M>E (q) 

xtr ^hc 


s--^ r 

rit 







FIG. 4: Qualitative behaviour of the metric function A(r) for 
the critical configurations (M = e{q)) in the cases Al (A(0) — > 
— CX3, case Al-4 in the text), A2a (A(0) = 1 — 16naq > 0, case 
A2a-2 in the text), A2b (A(0) = 1 - 16-Kaq < 0, case A2b-4 
in the text) and A2c (A(0) = 1 — 16naq = 0, case A2c-1 in 
the text). The three curves associated to each (A2) case cor- 
respond to the three ranges of the parameter a = 1, which 
determine their slopes at r = (see Eq. (|47|l l. The small frame 
displays the metric function for different values of the ADM 
mass in the case A2b with a > 1. The dashed line in this 
frame corresponds to the critical configuration (M — e{q), 
case A2b-4 in the text). The upper curves (M < e{q)) cor- 
respond to naked singularities (A), extreme black holes (B) 
and two-horizon black holes (C and D). The lower curves (E 
and F, with M > e{q)) correspond to the single- horizon black 
holes. Curves D and E show the approximation to the crit- 
ical configuration (for r > 0) in the limits M — e{q) — > 0^, 
respectively, illustrating the formation of the isolated discon- 
tinuities of the limit configurations at r = 0. The horizon 
radii for the extreme (rhextr) and critical (rhcrit) black holes 
are indicated. Similar figures to the one of the small frame 
can be straightforwardly displayed for the other values of o 
and for the cases Al, A2a and A2c, illustrating the discussion 
of section IV. 

As a well known example of this family let us consider 
the case of the original Born-Infeld model (whose gravi- 
tating ESS solutions were studied in Ref.jSi]), defined by 
the lagrangian density 



Lbi^Vbi{X,Y) 



(48) 



where is a constant, assumed positive. As — 
this function reduces to the Maxwell lagrangian density 
(Pm{X) = X. It is straightforward to check that this 
model satisfies the admissibility conditions of section II. 
The ESS solutions take the form 



Eir,q) 



(49) 



which are asymptotically coulombian, reduce to the 
Coulomb field as /i — ^ and behave around the cen- 
ter as E{r,q) ^ ^ — 27?^ ■ Thus the parameters of the 
expansion (1111) are a = l//i (the maximum strength of 
the field), ct = 4 and b = {2^^q^)~^ (which satisfies the 
relation (HH)). The energy density for this solution reads 



(50) 



and the exterior integral of energy is given by 



V 



1 1 5 

4' 2' 4' 



? 2 

-Mr 



(51) 



where Fi {x, y, z, u) is the Gaussian hypergeometric func- 
tion. As a consequence of the admissibility conditions, 
SexiT, 9) is a monotonically decreasing and concave func- 
tion for fixed q, which can be directly verified from its 
explicit expression. The total flat-space energy of the 
field takes the form 



(52) 



where T{x) is the Euler integral of second kind. 

If we solve now the equation (|23p for the extreme black 
holes we are lead to 



rhextr (q) 



(167rg — fi){Wnq + ^) 



327r 



(53) 



which requires the condition IQirq > /i to be fulfilled, 
in agreement with the preceding analysis (solutions A2b 
and A2c). The masses of these extreme black holes can 
be now immediately found by using Egs. ip^ . ((5T|) and 
([53| . As stated for general A2b solutions, the radii and 
masses of these black holes are increasing functions of 
the charge. The radii range from to c» as g ranges 
from TF- to cxD, whereas the masses range from ^(t^) 

167r ' ^ ^ 167r 

(extreme black point solution) to 00 in the same range of 

The verification of the remaining general statements 
for this BI example is now straightforward. 
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V. CONCLUSIONS AND PERSPECTIVES 

Let us summarize the main results found until now. 
We have considered the full family of NEDs whose la- 
grangian densities are functions of the two gauge invari- 
ants, restricted by some admissibility conditions and ex- 
hibiting asymptotically vanishing ESS solutions. When 
minimally coupled to gravity, we have established that 
the central field behaviour of these solutions (or, equiv- 
alently, the behaviour of the lagrangian densities on the 
boundary of their domain of definition at y = 0) fully 
characterizes the structure of their gravitating versions. 
In the case of NEDs with energy-divergent ESS solu- 
tions in flat space, the behaviour of the G-ESS solu- 
tions is qualitatively the same as the one of the RN solu- 
tions of Einstein-Maxwell field equations. In the case of 
NEDs with (flat-space) finite-energy ESS solutions, qual- 
itatively different features appear. Aside from naked sin- 
gularities, extreme black holes with a single horizon and 
two-horizons black hole solutions (as in the RN case), 
there are now non-extreme single-horizon black hole so- 
lutions and black-point solutions. In all cases there are, 
at most, two horizons. The presence of the "critical" so- 
lutions must be stressed. They can lead to naked singu- 
larities, as in the case (A2a-2), single- horizon black holes, 
as in cases (Al-4) and (A2b-4) or extreme black points 
as in the case (A2c-1). In all cases, the sequences of so- 
lutions (depending on the two integration constants M 
and q) converge, as M — > £((7), at fixed q, to limit metrics 
which coincide with the critical solutions for r > 0. Nev- 
ertheless, their structures at the center differ by singular 
((5-like or step) terms. This behaviour, which is similar 
to the one encountered in the Schwarzschild and RN so- 
lutions (23| . deserves a more careful analysis in order to 
elucidate the nature of the eventual point-like sources of 
the diverse solutions. 

At the center of these solutions there is always a curva- 
ture singularity, as established by the Bronnikov theorem 
[3 [i3j HO] • In fact this theorem applies only for asymp- 
totically coulombian fields (case B2). Nevertheless it can 
be easily generalized to the case Bl. The argumenta- 
tion of the theorem, as presented in Ref . [20j . fails in the 
case B3, but the regularity of the manifold at the cen- 
ter of the solutions can be excluded also in this case, for 
admissible models, by an alternative argument involving 
the positivity of the energy (see Ref. [231 ) • 

Asymptotically, if the fields behave as E{r) ~ (with 
p < —1, in order to endorse the convergence of the in- 
tegral of energy Sexir, q)), the corresponding metrics be- 
have as 



and approach asymptotic flatness as the Schwarzschild 
metric, excepting for solutions with vanishing ADM mass 
(which are always naked singularities) for which the grav- 
itational field (1 — A(r)) vanishes asymptotically at the 
same cadence as E{r). Obviously, solutions with p ~ —2 
(case B2) exhibit the same asymptotic behaviour as the 
RN one. 

Beyond the RN solution and the gravitating EH and 
BI models, employed here as representative of the fam- 
ilies of admissible theories in this context, the elabora- 
tion of other explicit examples reduces to the finding of 
the lagrangian densities (p{X, Y) satisfying the conditions 
displayed in Fig.l. This method, developed in Ref. [2l| . 
exhausts the class of this kind of admissible models. 

The thermodynamic analysis of Schwarzschild and RN 
black holes can be generalized to the present models. The 
results of Ref. [1^ show how the laws of thermodynamics 
can be consistently formulated for non-linear electrody- 
namic black holes. A discussion of the thermodynamic 
features of the gravitating NED models considered here is 
beyond the scope of this paper. This issue will be tackled 
elsewhere (25| . 

Finally, some comments on the stability of the G-ESS 
solutions are in order. The necessary and sufficient condi- 
tion for the linear stability of the ESS solutions of NEDs 
in flat space has been obtained in Ref. 21| and requires 
the inequality 



dX- dY^' 



(55) 



A(r) 



2M 



32Trq 



(p-2)(p+l) 



(54) 



to be satisfied by the lagrangian densities in the domain 
of definition of the solution. It is not clear a priory 
whether the effect of the gravitational field improves 
the stability of solutions which are unstable in fiat 
space, or whether new conditions should be added to the 
lagrangian densities of the admissible NEDs supporting 
stable ESS solutions in fiat space, to maintain stability 
of their gravitating generalizations. The analysis of this 
issue is in progress j27| . 
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